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Abstract. The single-particle spectrum of deformed shell-model states in nuclei, is shown 
to exhibit a supersymmetric pattern. The latter involves deformed pseudospin doublets and 
intruder levels. The underlying supersymmetry is associated with the relativistic pseudospin 
symmetry of the nuclear mean-field Dirac Hamiltonian with scalar and vector potentials. 



1. Introduction 

The concept of pseudospin doublets [TJ [2] is based on the empirical observation in nuclei of quasi- 
degenerate pairs of certain normal-parity shell-model orbitals with non-relativistic single-nucleon 
quantum numbers 

(n,i,j = £ + 1/2) and (n - 1,£ + 2, j = £ + 3/2) . (1) 

The doublet structure is expressed in terms of a "pseudo" orbital angular momentum, 1 = 1 + 1, 
and "pseudo" spin, s = 1/2, which are coupled to j = £ + s. For example, (W5/2, 0g 7 / 2 ) and 
(2si/ 2 ) 1^3/2) wm have I = 3 and I = 1, respectively. The states (n = 0,£,j = £ + 1/2), 
with aligned spin and no nodes, are not part of a doublet. For large j they are the "intruder" 
abnormal-parity states, i.e., 0g 9 / 2 , 0/i n / 2 , 0i 13/ / 2 , which are unique in the major shell. In the 
presence of axially-symmetric deformation, the doublets persist [3] with asymptotic (Nilsson) 
quantum numbers 

[N,n 3 ,Ap = A + 1/2 and [JV, ra 3 , A + 2]0 = A + 3/2 . (2) 

Here N, 713, are harmonic oscillator quantum numbers, A and are, respectively, the components 
of the orbital and total angular momentum along the symmetry (3rd) axis. In this case, the 
doublets can be expressed in terms of pseudo-orbital, A = A + 1, and pseudospin, ft = +1/2, 
projections, to obtain fl = A+l/2. For example, ([420] 1/2, [422]3/2) and ([411]3/2, [413]5/2) will 
have A = 1 and A = 2, respectively. The intruder nodeless states, [N, n^, A = N— n^]Q = A+l/2, 
are not part of a doublet. Pseudospin doublets play a central role in explaining features of 
nuclei [3J , including superdeformation [4J and identical bands [5j [6] . They have been shown to 
originate from a relativistic symmetry of the Dirac Hamiltonian in which the sum of the scalar 
and vector mean-field potentials cancel [7]. The generators of this relativistic SU(2) symmetry 
have been identified [8] and consequences for radial nodes have been considered [9]. Recently, it 
has been shown that a pseudospin-invariant Dirac Hamiltonian with spherical potentials, gives 
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Figure 1. Nuclear single-particle spectrum 
composed of deformed pseudospin doublets and an 
intruder level. All states share a common A. 
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Figure 2. Typical spectra of good 
and broken SUSY. The operators 
and Q + connect degenerate 
states in the i2+ and H- sectors. 



rise to a supersymmetric pattern [TO] . In the present contribution we extend this study and show 
that a supersymmetric pattern occurs also when the relevant potentials are axially deformed [11] . 

Fig. Q] portrays the level scheme of a subset of deformed pseudospin doublets, Eq. ([2]), with 
fixed N, A, f2, f2' and n± = (N — A — n^)/2 = 1,2,3,... together with the intruder level 
[N,ri3, A = N — ri3]Q, = A + l/2, (or equivalently, n± = 0). The single-particle spectrum exhibits 
towers of pair-wise degenerate states, sharing a common A, and an additional non-degenerate 
nodeless intruder state at the bottom of the spin-projection aligned tower. A comparison with 
Fig. 2 reveals a striking similarity with a supersymmetric pattern. In what follows we identify 
the underlying supersymmetric structure associated with a Dirac Hamiltonian possessing a 
relativistic pseudospin symmetry. 



2. Dirac Hamiltonian with axially-deformed potentials 

A relativistic mean-field description of nuclei employs a Dirac Hamiltonian, H = a • p + 
/3(M + Vg) + Vy for a nucleon of mass M moving in external scalar, Vs, and vector, Vy, 
potentials. When the potentials are axially-symmetric, i.e., independent of the azimuthal 
angle (j), Vsy = Vsy(p,z), p = \/x 2 + y 2 , then the z-component of the angular momentum 
operator, J z , commutes with H and its half-integer eigenvalues Q, are used to label the Dirac 
wave functions 



( g+(p,z)e^- l l 2 ^ \ 
g-{p,z)e^ +1 l 2 )<t> 
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(3) 



Here = g ± (p,z) and / ± = / ± (p, z) are the radial wave functions of the upper and lower 
components, respectively. The Dirac equation, = E^, leads to a set of four coupled partial 
differential equations for these components. For each solution with Q > 0, there is a degenerate 
time-reversed solution with — £1 < 0, hence, we confine the discussion to solutions with 0, > 0. 
Of particular interest are bound Dirac valence states with < E < M and normalizable wave 
functions. The potentials enter the Dirac equation through the combinations 



A(p,z) = E + M + V s (p,z)-Vy(p,z) 
B(p,z) = E-M-V s (p,z)-Vy(p,z) 



(4) 



For relativistic mean-fields relevant to nuclei, Vs is attractive and Vy is repulsive with typical 
values Vs(0) ~ -400, Vy(0) ~ 350, MeV. The potentials satisfy pV s (p, z), pVy(p, z) -> for 



p — > and Vs{pi z), Vv(p, z) — > for p — > oo or z — > ±00. The boundary conditions imply that 
the radial wave functions fall off exponentially for large distances and behave as a power law for 
p — s> 0. Furthermore, for z = and p ->• 00, f~/g + oc (M - E) > and g~ / f + oc (M + > 0, 
while for z = and p — >■ 0, f~/g + oc B(0)p and g~ / f + oc —^4(0)/). These properties have 
important implications for the structure of radial nodes. In particular, it follows that for 
potentials with the indicated asymptotic behaviour and A(0), B(0) > 0, as encountered in 
nuclei, a necessary condition for a nodeless bound eigenstate of a Dirac Hamiltonian is [11] 

g - =0 or /+ = . (5) 



3. Relativistic pseudospin symmetry in nuclei 

A relativistic pseudospin symmetry occurs when the sum of the scalar and vector potentials is 
a constant [7j 



V s (p,z) + V v (p,z) = A . 



(6) 



The symmetry generators, Si, commute with the Dirac Hamiltonian and span an SU(2) 
algebra [5] 
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Here Sj = o"i/2 are the usual spin operators, defined in terms of Pauli matrices. The Dirac 
eigenfunctions in the pseudospin limit satisfy 
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(8) 



and form degenerate SU{2) doublets. Their wave functions have been shown to be of the 
form [T2] 
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(9) 



where A = f2 — /I > is the eigenvalue of J z — S z . Near the pseudospin limit, eigenstates with 
A ^ occur in degenerate doublets. An exception to this rule are eigenstates whose radial 
components have no nodes and, therefore, are not part of a doublet. The latter property follows 
from the fact that a nodeless bound Dirac state satisfies the criteria of Eq. ([5]), hence has a wave 

function V?^ _jf_i/2' as m ® ' w ^ n 5 + > ff> / 7^ an d //<7 + > 0. Its pseudospin partner state 
has a wave function *1 _l ,„• The radial components satisfy Bg~ = [B — 2(A / p)f /g + ]g + , 

where -B is defined in Eq. This relation is satisfied, to a good approximation, for mean- field 
potentials relevant to nuclei, and the r.h.s. is non-zero and, consequently, g~ ^ 0. If so, then the 
partner state, -A+1/2' * s a ^ so n °deless, but it cannot be a bound eigenstate since its radial 
components do not fulfill the condition of Eq. ©. 

The relativistic pseudospin symmetry has been tested in numerous realistic mean-field 
calculations of nuclei and was found to be obeyed to a good approximation, especially for 
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Figure 3. Grouping of deformed shell-model states [N = 4, n3,A]f2, exhibiting a pattern of 
good SUSY, relevant to the pseudospin symmetry limit [TT] . 

doublets near the Fermi surface [121 021 [IH EH IS [TTJ, [18]. The dominant upper components 
of the states in Eq. ([9]), involving g + and g~, correspond to the non-relativistic pseudospin 
doublets of Eq. ([2|). The nodeless Dirac eigenstates correspond to the non-relativistic intruder 
states, [N,ri3,A = N — 713]^ = A + 1/2, mentioned above. In what follows we show that the 
ensemble of deformed pseudospin doublets and intruder levels, forms a super symmetric pattern 
whose origin can be traced to the relativistic pseudospin symmetry of the Dirac Hamiltonian. 



4. Relativistic pseudospin symmetry and supersymmetric quantum mechanics 

The essential ingredients of supersymmetric quantum mechanics [TU] are the supersymmetric 
Hamiltonian, T~L, and charges Q+, Q- = Q\, which generate the supersymmetry (SUSY) algebra 

[H,Q±] = , {Q±,Q±} = , {Q-,Q+} = n. (10) 
Accompanying this set is a Hermitian .^-grading operator satisfying 

v ] = v , [h,v] = o , {Q±,v} = o , v 2 = t. (11) 

The +1 and —1 eigenspaces of V define the "positive-parity", H + , and "negative-parity", H-, 
sectors of the spectrum, with eigenvectors vpW and respectively. The SUSY algebra 

imply that if ^(+) is an eigenstate of H, then also = [ s an eigenstate of % with the 

same energy, unless Q-Vl'W vanishes or produces an unphysical state, (e.g., non-normalizable) . 
The resulting spectrum consists of pairwise degenerate levels with a non-degenerate single state 
(the ground state) in one sector when the supersymmetry is exact. If all states are pairwise 
degenerate, the supersymmetry is said to be broken. Typical spectra for good and broken SUSY 
are shown in Fig. [2j 

Degenerate doublets, signaling a supersymmetric structure, can emerge in a quantum system 
with a Hamiltonian H, from the existence of two Hermitian, conserved and anticommuting 
operators, R and B 

[H,R] = [H,B] = {R,B} = . (12) 

The operator R has non-zero eigenvalues, r, which come in pairs of opposite signs. B 2 = 
WB = f(H), is a function of the Hamiltonian. A ^-grading operator, V r , and Hermitian 
supersymmetric charges, Q\, Q2 and Hamiltonian H can now be constructed as 

V r = R/\r\ , Qx = B , Q 2 = iQiPr , U = Q\ = f(H) . (13) 

The triad of operators Q± = (Qi ± iQ2)j2 and % form the standard SUSY algebra of Eq. (TlOj) , 
Focusing the discussion to the relativistic pseudospin limit, the elements of the SUSY algebra 

are 

Q 1 = 2(M + A -H)S X , Q 2 = 2(M + A -H)§ y , V r = 2S z , 

H = Qj = (M + A -H) 2 , (14) 



Table 1. Conserved, anticommuting operators, R and B, for Dirac Hamiltonians (H), with 
axially-deformed scalar and vector potentials, exhibiting a super symmetric structure [11]. Here 

Si = ^ u p s ^ u p Q ^J and 5i = ( [/ s° [/ ) are ^ e P seu dospin and spin generators, respectively. 
R z = [ M + V s (z) ] ± 3 + 15P Z and R p = [ M + Vs(p) ]E 3 -^ 75 (Sxp) 3l where = (<£ ° ) . 
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Table 2. As in Table Q] but for spherical scalar and vector potentials [101 [20] in the Dirac 
Hamiltonian (H). Here K = -0 ( £ • £+ l), with £ = r x p, f z = (S x L) z , f z = (S x 

y = L • L + 1/4 — jj and Y = L • L + 1/4 — j|. The operators Sj, Sj and Si are defined in the 
caption of Tableland L = Ji — Si, L = Ji — Si. 



Dirac Hamiltonian 
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where H is the pseudospin-invariant Dirac Hamiltonian with potentials satisfying the condition 
of Eq. (|6j). The operator B = Q\ connects the doublet states of Eq. Q. The spectrum, for 
each A / 0, consists of twin towers of pairwise degenerate pseudospin doublet states, with 
£l\ = A — 1/2 and CI2 = A + 1/2, and an additional non-degenerate nodeless state at the bottom 
of the Q\ = A — 1/2 tower. Altogether, the ensemble of Dirac states with — ^1 = 1 exhibits 
a supersymmetric pattern of good SUSY, as illustrated in Fig. [31 

The supersymmetry discussed above, is one of several possible classes of supersymmetries 
of Dirac Hamiltonians in 3+1 dimensions, with cylindrically-deformed scalar and vector 
potentials [11]. The supersymmetries associated with the pseudospin and spin limits arise 
when the potentials obey a constraint on their sum or difference, respectively. Additional 
supersymmetries arise when one of the potentials has a 1/p dependence and the second potential 
depends on z. The relevant R and B operators (|12j) of the various supersymmetries, are listed in 
Table[T] A similar analysis can be done for spherical scalar and vector potentials and, as shown in 
Tabled the corresponding supersymmetries are associated with the pseudospin (Vs+Vy = Aq), 
spin (Vs — Vy = Sq) and Coulomb (Vs, Vy cx 1/r) limits of the Dirac Hamiltonian [TOl [20], 



5. Summary 

As previously noted, the single-particle spectrum of spherical shell-model states in nuclei, 
exhibits characteristic features of a super symmetric pattern. The latter consists of twin towers of 
pair- wise degenerate spherical pseudospin doublets (with n = 1,2...) sharing a common £, and 
an additional non-degenerate nodeless (n = 0) intruder state at the bottom of the spin-aligned 
tower 

(n, I, j = I + 1/2 ) (n - 1, £ + 2, j = £ + 3/2 ) 

(n = 0, 4 j = £ + l/2) 

3 = 1-1/2 j = £ + l/2 

As discussed in the present contribution, a similar situation is observed in the single-particle 
spectrum of deformed shell model states. The corresponding supersymmetric pattern consists of 
twin towers of pair- wise degenerate deformed pseudospin doublets (with nj_ = 1,2...) sharing 
a common A, and an additional non-degenerate nodeless (n± = 0) intruder state at the bottom 
of the spin-projection aligned tower 

[N, n 3 - 2nj_, A] = A + 1/2 [N,n 3 - 2n ± , A + 2] Q = A + 1/2 
[N,n 3 ,A]Q = A + 1/2 

ft = A- 1/2 ft = A + l/2 

with A = N — n 3 . In both cases, the underlying supersymmetric algebra is associated with the 
relativistic pseudospin symmetry of the nuclear mean-field Dirac Hamiltonian, involving scalar 
and vector potentials. Additional supersymmetries of Dirac Hamiltonians with such mixed 
Lorentz structure can be identified. 



n = l,2,... 

n = 



(15) 



n ± = 1,2,.. 
n ± = 

A = A + l 



(16) 
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